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A theory for an electron affinity of ionic clusters is proposed both in a quasiclassical approach and 
with quantization of a polarization electric field in a nanoparticle. An interaction of an electron with 
longitudinal optical phonons in an ionic cluster is described. A critical size of the cluster regarding in 
formation of an electron's autolocalized state, dependencies of an energy and a radius of a polaron on 
a cluster's size are obtained by a variational method. It has been found that a binding energy of the 
electron in a cluster depends on a cluster's radius but the radius of electron's autolocalization does 
^ , not depend on cluster's radius and equals to the polaron radius in a corresponding infinity crystal, 

Q ■ moreover the bound state of the electron in a cluster is possible if the cluster's radius is more then 

the polaron radius. A perturbation method for finding of the critical parameters is proposed, where 
an dependence of the critical size on a momentum of an electron in a conduction band is observed. 
. Keywords: cluster, polaron, autolocalization, phonon confinement. 
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C/2 ! I. INTRODUCTION 

<D 

A bound state of an electron in a polar media is caused by a local polarization of the medium inducted by the 
electron. The electron is in a potential well in states with a discrete energy and its field supports the polarization of a 
lattice. An interaction of electrons with optical phonons resulting in an instability of an electron's state in a conduction 
band is described by both a quasiclassical approach and an electron-phonon Frohlich hamiltonian. On basis of belief 
about a self-consistent state of an electron and a polar medium (ammonium, water, protein globules) mathematics 
models of an excess electron in a cluster have been built and investigated [U-Q . In the models it is supposed a solvated 
electron does not belong to a separate molecule but collectively interacts with many atoms of the polar medium. It is 
O ■ shown by analytical and numerical methods that a critical size of molecular clusters consisting of polar molecules (for 
example (H20)„ or (NH3)„) exists. In a cluster with a less size than the critical size the autolocalized state (polaron 
I . state) can not exist. 

" Since the polaron effect exists in an ionic crystal then we can suppose the effect can occur in an ionic cluster too 
(for example Nai4Cli3) beginning with a some critical size. In the papers [7H9| electron localization in alkali-halide 
clusters M+X~_ x and M+X~_ 2 was numerically investigated. It was supposed that the clusters contain F-center 
defects and electron localized near an halide ion. Obviously, in this model an effect of localization of an electron in 
a cluster takes place, but not autolocalization because main contribution in a binding energy of the excess electron 
t-H 1 is Coulomb-like interaction with F-center. The critical size in this case is not observed. We will consider the system 
neutral ionic cluster+excess electron where we have interaction of the electron with an induced polarization only and 
process of autolocalization. 

The system electron+polarization can be described in a quasiclassical approach. This means a motion of an electron 
is quantized (is described by Schrodinger equation), however specific polarization of a medium P(r) is considered as a 
classical variable and described by electrostatic equations. However an interaction of a electron with the polarization 
field can be described in quantum approach. Borders of a cluster essentially influence upon quantization of optical 
oscillations - the optical phonon confinement [lfj , that is not considered in the quasiclassical approach. Influence of a 
geometric upon q uantization of optical oscillations and an interaction of an electron with them have been investigated 
in papers |lll4l3| | for quantum wells, wires and boxes. However, it should be noted that polaron states in the 
heterostructures essentially depend on a confining potential. But in a nanoparticle the polaron effect is stipulated 
by interaction with an induced polarization only. A basic difference of optical and acoustical modes in a cluster and 
heterostructure from the modes in an infinity crystal is a restriction to possible wavelengthes of the oscillation and 
presence of interfacial phonons. These properties of phonons cause differences of a polaron state in a nanoparticle 
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from the polaron state in an infinity crystal. 

In a section [TT] wc determine an energy of a polaron in a nanoparticle and critical parameters of nonoparticles 
regarding to formation of a polaron state by a variational method in a quasiclassical approach. The variational 
method, unlike numerical calculations, lets obtain analytical expressions for connection of an autolocalized state's 
energy and a polaron radius with parameters of a cluster (a size, an electron phonon coupling constant, a frequency of 
lattice optical oscillations). In a section [TTT] a field of deformations of an ionic cluster is quantized taking into account 
the boundary condition on its surface. The energy of a polaron state, the critical size of a nanoparticle and the 
polaron radius are determined by a variational method on the basis of an electron-phonon hamiltonian for a spatially 
limited medium. Obtained results is compared with result obtained by the numerical methods earlier. In a section 
IIVI a perturbation theory method was proposed in order to find the critical parameters of a cluster. 



II. THE QUASICLASSICAL APPROACH. 

In this section we will consider a quasiclassical description of a electrons' interaction with a polarization displacement 
of ions from their equilibrium positions. Let us suppose an electron is localized in an ionic cluster of a spherical shape 
with a radius R and described by a wave function ^(r). A created by electron mean Coulomb field divD(r) = e|\l/(r)| 2 
induces a polarization of the cluster. In turn the polarization electric field acts on the electron. A dipole moment of 
cluster's unit volume P(r) is determined by a difference of static and high-frequency polarizations: 

P(r)=P (r)-P oo (r)= (— --\t>=Ix>, (1) 

\£oo ej e 

where D is an electric displacement, e and are static and high-frequency permittivities accordingly. For a medium 
surrounding the cluster (for example the vacuum) we have — e. Then 

1 = £^ r<R 

(2) 

e = 00 r > R 

A polarization of a cluster brings to appearance of a density of polarization charge p = — divP. Let E = — \7tp is 
an electric field of the polarization, moreover divE = p/ea =>■ D = — £o£E, where £0 is an electric constant. Then 
an energy of a cluster's polarization is [TEl ]: 

Afield = \J = \e* J ^ 2 dV. (3) 

Hence an energy functional of a system electron+polarized cluster has a form: 

<P) = ^J \W\ 2 dV + e J \*\ 2 <pdV + i £o J e(V<pfdV, (4) 

where the first term is a kinetic energy of a localized electron in the state Vf(r), the second term is an interaction 
energy of an electron with the electric field ip(r) of the induced polarization. A localized state of an electron in a 
cluster is energetically profitably if / < 0. Varying the functional (0| with respect to "J at a condition J \^\ 2 dV = 1 
we obtain an equation describing movement of an electron with an energy E in a potential well etp created by the 
electron: 

h 2 

A* + eip^J = EV. (5) 

2m 

Varying the functional (01 with respect to (p taking into account © we obtain an equation describing a field of the 
polarization induced by the electron: 



A^^si^t r<R 
A(p = r>R 

Boundary conditions on a surface of the cluster (1) and the environment (2) have a view: 



(6) 



Vi = V2, ^=^, (7) 
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where n is a normal to the cluster's surface. Since £2 = 00 and in view of Eq. 
of the cluster we can suppose 



ip = 



for 



r> R. 



I, on a border of the cluster and out 

(8) 



The boundary condition ([SJ tp{R) — is a main difference of an electron's autolocalization in a cluster from a polaron 
state in a infinity crystal where the electrical field of a polarization is equal to zero at the infinity only: </?(oo) = 0. 

To obtain energies of polaron states and a critical size of a cluster we will use a variational method because it lets 
find analytical expressions. Let the cluster is characterized by values of static e and high-frequency £00 permittivities 
and effective mass to of an electron in a crystal's conduction band (for NaCl to = 2.78mo, too is a electron's mass in 
vacuum). Obviously these values have sense for clusters with sizes which are much more than interatomic distances 
(~ 2 -j- 3A) only. However a proposed below method can be extrapolated to small clusters. With help of the equations 
© and the condition <p{oo) = the functional (0} can be simplified: 



2m 



[\W\ 2 dV-^s Jl(V<p) 2 dV. 



(9) 



For a region out of the cluster r > R we suppose £(p — 0. A wave function of a system's (electron+polarized cluster) 
ground state can be selected in a form [la ]: 



* = 



1 + r /r 



■exp 

ro 



(10) 



where the radius ro is a variational parameter - a polaron radius. As far as the cluster is spherical, then the equation 
© can be rewritten in a view: 



r 2 dr 
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The first integral and a solution of this equation are functions: 
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y\ 2 r 2 dr + d 

m 2 r 2 dr + d ) <lr + C-_ 



(11) 

(12) 
(13) 



where the constant C\ must be selected for ip(0) 7^ 00, the constant Ci must be found from the boundary condition 
©. Substituting the wave function (fTU|) in Eqs. (ll2ll3[) we have: 
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Substituting ^ (ITU1) and ^ (|14ll6p in the functional /(vf, <p) © we obtain an energy of a polarized cluster with a 
localized electron as a function of a polaron radius: 
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The following terms give a main contribution in the energy at tqR > 0.5: 
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(18) 



where ioo is an energy of an infinity crystal with a localized electron. An autolocalized state of an electron in a cluster 
(or in a crystal) is energetically profitable if 7 < 0. From Eqs. (|17ll8|) we can see a spatial limitation of a crystal 
increases the energy of the polaron state: 7^ < I. From Eq. (|18p we can see the polaron's radius r does not depend 
on the cluster's radius R: 
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(19) 



The function J(fo) (H3 f° r different values of the cluster's radii R is shown in Fig. ((TJ) . We can see, that a bound 
state of an electron in a cluster consisting of molecules NaCL can exist beginning from R sa 3a, where a = 2.81A is an 
interatomic distance. With increase a cluster's size a bound energy increases and aspires to a polaron energy in an 
infinity crystal loo ~ — 0.265eV. A radius of electron's autolocalization ro in a cluster minimizes the energy I. From 
Fig. (JTJ) we can see ro does not depend on a cluster's radius almost and for sodium chloride the radius is equal to 
ro = 1.5A. The dependencies shown in Fig.([T]) is well approximated by the more simple dependence (|18p . 
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Figure 1: The dependencies of an energy of a system electron+cluster (|17|) on a polaron radius I = I(tq) for an ionic cluster 
consisting of molecules NaCL. The curves have been built for cluster's radii: 2a, 3a, 8a, 15a, 40a, where a = 2.81A is an interatomic 
distance. The bold line is the energy of an infinity crystal with an electron. 

The equations (I5I6|) were solved numerically in a paper Q for clusters of ammonia and water. It was shown that 
critical size of the clusters consisting of polar medium exists. However analytical method discovers the next pattern. 
The polaron radius in a cluster does not depend on a radius of the cluster. The bound state of an electron in the 
cluster exists till a polaron radius is less than the cluster's radius. That is the autolocalized electron is hidden in 
depth of the cluster always and has a negligibly low probability to be on the outside. If the cluster is enough small so 
that the electron can be out of the cluster ,the bound state disappears. Above-mentioned particulars are illustrated 
in Fig.@. Similar result ware obtained in [J, but another functional I^, f) was considered and ground state energy 
of an electron in a cluster and the autolocalization radius were expressed by static dielectric constant only. Hence in 
this model autolocalisation effect is not polaron effect unlike our model. 

Since a polaron radius ro is independent on a cluster's radius R, then we can obtain an energy of a system 
electron+cluster as a function of the cluster's radius I(R) assuming in the formula (|17[) the polaron radius is equal to 
the optimal value (fT9"]). The result is shown in Fig. ([3]), where we can see the critical radius of a cluster (consisting of 
NaCL) determined by the equation I(Rcr) = is equal to i? C r = 7. OA (2-^3 interatomic distances). With increase of 
the cluster's radius the energy asymptotically aspires to its value in an infinite crystal I(R —¥ oo) — > Joo- 

A ground state of an electron in the potential well eip is Is. The energy /('J, <p) (|4l9p is a sum of a ground state 
energy of an electron in the potential well and an energy of a cluster's polarization l/2eo J e(V(f) 2 dV > 0. Then the 
ground state energy of an electron in a cluster is 

E(R) = J |W| 2 dV -e J I(Vf) 2 dV, (20) 

moreover the wave function "J and the field if is taken in the forms (|10l) and (|14M16I) accordingly. The radius r 
is equal to its optimal value (|19[) . The energy E(r) is a necessary energy to quickly remove an electron from the 
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Figure 2: The probability distribution of being of an electron in a cluster. The vertical lines means radii of clusters for which 
an energy of an electron's autolocalized state was calculated in Fig.fT). The radius 2a is subcritical (the bound state is absent), 
the radius 3a is some bigger than the critical radius - the bound state appears. 
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Figure 3: The energy of a system electron+cluster (consisting of NaCL) as a function of a cluster's radius I(R) (the solid line). 
The doted line is an energy of a polaron in an infinite crystal R = oo. Rcr is the critical radius of a cluster. 



localized state to a free state (for example by absorption of a photon with an energy frui > \E\). The graph E(R) 
is shown in Fig.Q, where we can see that ls-state is energetically profitable beginning from a radius of the cluster 
R S CT sa 4A < Rcr = 7 A (substance NaCl). That is in the interval R CT < R < Rcr a metastable autolocalized state of 
an electron in an ionic cluster is possible. 

If the potential well eip (fT5|) is sufficiently deep then other discrete energy levels are possible. At increasing of the 
cluster's radius p-state can appear, s — p photo-transitions occur without changes of ions' positions - FranckCondon 
principle. That is the electrical field ip of a polarization and a polaron radius ro are the same as the values in s-state: 
(fT5| and (fT§|) accordingly, however a wave function of an electron can be chosen in a form: 

* = f— V — exp f-— Vosfl, (21) 
\r J nr \ r J 

where £ is a variational parameter, 9 is a polar angle. Then an energy of an electron in p-state is: 

E p = - / 2irr 2 dr sin 6d9 |W p (r, 6)\ 2 + e / 2nr 2 dr sin BdB |* p (r, 9)\ 2 <p{r). (22) 

2m* J Jo J q J 

The expression (|22]l is a function of a variation parameter £. The graphics of a dependence E p on £ for clusters' 
radii R = 2a, 3a, 8a, oo are shown in Fig.([S]), where we can see an optimal value (it minimizes E p ) of the parameter £ 
weakly depends on a cluster's radius and it can be supposed £ = 0.65. Then we can easily find a dependence of p-state 
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Figure 4: The energies of a electron's ground state Is and an electron's p-state in the potential well eip as functions of a cluster's 
radius E(R) and E P (R) (solid lines). The dotted lines show the energies of an electron in an infinity crystal R = oo. Rct is 
a critical radius when the polaron state is energetically profitable (see Fig[3]), i?| r and R^ r are a critical radii when ls-state 
and p-state of an electron in the potential well (|15I16[) is energetically profitable accordingly. 

energy on a cluster's radius E p (R) - Fig.(j4]), where we can see p-state is energetically profitable beginning from the 
radius of a cluster R P CT « 8A > R C z > Rcr- With increasing of a cluster's radius the energy E p (R) asymptotically 
aspires to its value in an infinity crystal. In a bigger clusters d, f and so on states are possible. 




Figure 5: The dependence of electron's energy in p-state in a cluster on a parameter £ for an ionic cluster consisting of NaCL. 
The curves have been built for radii of clusters: 2a, 3a, 8a, where a = 2.81A is an interatomic distance. The bold line is electron's 
energy in p-state in an infinite crystal. 

It necessary to notice that all found critical sizes for the cluster have values of the order of several interatomic 
distances a. This means the continuous medium model is not good for description of a such nanoparticle. Moreover 
an effective mass to in a conduction band figures in the model. Hence the model is correct for enough big clusters 
R a. Sodium chloride is a strongly polarizable substance. For less polarizable substances the model can be used 
with a better accuracy. 

An essential difference of a polaron state of an electron in a cluster from the polaron state in an infinity crystal is 
localization of an electron in the cluster's center in result of the boundary condition (|5|). On the contrary in an infinity 
crystal all points are equivalent (in a continuous representation of a substance) therefore the polaron can move. In 
this case the polaron is characterized by an effective mass M 3> to. 
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III. THE QUANTUM APPROACH. A STRONG COUPLING THEORY OF AN ELECTRON IN A 

CLUSTER. 



An interaction of a electron with the polarization field can be described in quantum approach. Borders of a cluster 
essentially influence upon quantization of optical oscillations that is not considered in the quasiclassical approach. A 
basic difference of optical and acoustical modes in a cluster from the modes in an infinity crystal is a restriction to 
possible wavelengthes of the oscillation. In a such situation it is necessary to quantize a field of a cluster's polarization 
tp. In a harmonic approximation a hamiltonian of a system electron+cluster has a form: 



^ A 1 ~ 

H= - — A + -e a e 
2m 2 



^(V^) 2 + (V^) 5 



dV + eip = T + U f + V it 



(23) 



where the first term T is an operator of a kinetic energy of an electron, Uf is an operator of a deformation energy 
of a crystal (cluster), Vi n t is &n operator of an interaction of the electron with the deformation field. The procedure 
of quantization of the deformation field lies in the fact that the electric field ip must be expressed in creation and 
annihilation b operators of optical phonons so that 



dV = hcoJ2 ( 6 q & q 



(24) 



where q is a wave vector of a phonon. For simplicity we will suppose optical oscillations is dispersionless u> = const. 
The field ip is proportional to a displacement of atoms £. The displacement and the field must satisfy periodic 
boundary conditions: 



£n — in \2i\- Pn — P>n+N2a, 



(25) 



where 2a is a lattice constant of a two-component ionic crystal. Then interaction of an electron with the deformation 
field (electron-phonon interaction) is written as follows: 



Vint = e& ■ 



Tioj 



x - e 



iqr 



2£n?V4^ 



(26) 



where q ia taken from the first Brillouin zone — 7r/2a < q x ,Qy,Qz < 7r/2a. The creation and annihilation operators 



are written in sense be 



fct 



9q — f oqe~~", Oq — f bqe lut . The cyclic boundary conditions are ensured by periodic multipliers 
e iqr_ p Qr an m g n j(;y crystal we can neglect discreteness of a substance and we can assume a period of a lattice is 
infinitely small, hence — oo < q x , q y , q z < oo. 

Due ionic type connection between atoms a shape of ionic cluster is close to cubic (for example unlike a ball 
shape of metallic cluster). Computer simulations and experimental evidences 04^1 speak about alkali halide clusters 
M+A"~, n ~ m > 14 form stable cuboid structures and have a cubic microlattice. Clusters as and nucleuses can have 
magic numbers of atoms - such numbers of atoms when all shells are filled. Magic clusters are most stable and magic 
ionic clusters have cuboid shape. For large cluster U>aa shape can any (a is an interatomic distance). Obviously 
for large clusters R ^> the shape have not principal importance because an excess electron has low probability to 
be near a surface of the cluster. For convenience we will consider the large cluster as cubic in this case too. 



Let we have an ionic cluster (for example consisting of NaCL) of a cubic form with sizes 



< x,y,z < 



A 



main difference of a cluster from an infinity crystal is the condition on a cluster's boundary (j5J) instead of the cyclic 
boundary condition (|25j) . that causes the phonon confinement - Fig 17] That is for the cube the boundary conditions 
have a view: 



ip I x = ±- 



L 



<p[y = ±—)=<p[z = ± 



= 0. 



(27) 



The field tp (|26|) does not satisfy the conditions (|27|) . In the article a method of construction of a hamiltonian for 
an electron-optical phonon interaction in a dielectric slab (with boundary conditions <p (z = ±-§) = 0) was proposed. 
But for a quantum box (cluster) the method brings to a very cumbersome expression. We propose the following 
method. The boundary condition is satisfied by the operator 



.4 



1v 



cos(q x x), odd 
sin(g x x), even 



cos(q y y), odd 
sin(aj,y), even 



cos(q 2 z), odd n z 
svn{q z z), even n z 



(28) 
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where projections of a wave vector q of optical phonon are q x 



L 

y-.z 2 



"x,y ,z 2 



(such that cos [q x 



y> z 2 / 



0, sin \Qx,y,zM) — 0). The unknown constant A is found from the condition of secondary quantization of the 

deformation field (IM1) . From where we obtain A = ../ -^U. Energy of an electron is minimal when it is situated in the 

center of a cluster and a probability amplitude is symmetrical about the center ^(r) = 4'(— r). Then the summands 
with sin may be omitted because they do not give any contribution in an electron-phonon interaction energy (and 



L 11 

hence in a polarization - the field ip): J_ L / 2 

dip 



^ (x) sin(q x x) = 0. This is equivalent to boundary conditions 

dip 



dip 



0) 



-(z = 0) = 0. 



(29) 



The field ip must satisfy the boundary conditions (|27I29I) on a cluster's surface. Whole number n of half-waves of 
phonons must be placed in a cube edge L — na. Moreover a half of phonon's wavelength must not be bigger than 
length of a cube edge and smaller than interatomic distance a (between a cation and an anion on a cube edge) 
a < A/2 = ~ < L. Hence the possible projections of wave vector q = {q x ,q y ,q z ) of an optical phonon must take on 
values: 
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Then a hamiltonian of a system electron+cluster is written in a view: 
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L 3 ^ q 



where 



M n 2 = e 2 



2e £ 
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e / m \ 



1/2 



(31) 



(32) 



Here a plays a role in a electron-phonon coupling constant. We suppose optical phonons are dispersionless and a 
does not depend on sizes of a cluster. Experimental data say about an oscillatory spectrum u)(q) for clusters (NaCL)jy 
quickly approaches with an oscillatory_spectrum of an infinity crystal at increase of N. So N — 4 is enough for having 
characteristics of an infinity crystal |14| . Since an autolocalized electron is hidden in a depth of the cluster always 
and has a negligibly low probability to be on the outside (it was shown in a previous section) then we can neglect of 
an interaction with interface phonons. 

Let us find specific polarization operator P using a set of equations: 



.£_ 

£0 



A<p 

p = — divP 

where p is a density of polarization charges. Then 



(33) 



huj 



2e eL 3 ^ q 

x [iq x sin(q x x) cos(q y y) cos(q z z) + jq y cos(q x x) sin(q y y) cos(q z z) + kq z cos(q x x) cos(q y y) sm(q z z)} . (34) 

In turn, the specific polarization vector is proportional to displacements of atoms of a crystal (cluster): P = v ^ . 
The fields ip and ^ can be qualitatively illustrated by the example of an atom chain. Then in ID system we have 
ip cx cos(qx), ^ cx q sin(qx). These configurations of the fields are shown schematically in Fig. ([?])• We can see in a chain 
consisting of four atoms (two anions and two cations) the oscillations are possible with two wave vectors q = ir/3a 
and q = it /a only in the course of the boundary conditions (|27I29[) . 

It is convenient to transit from the creation and annihilation operators to new conjugate operators: 



Qq 



V2 



Pa 



\/2 



ht 



(35) 



Figure 6: The schematic illustration of possible oscillations (in the course of the boundary conditions (|27l29[l ) of an electric 
polarization field <p and an atom displacement £ for a chain consisting of two cations and two anions . 



Then the hamiltonian is written in a view: 

« = -^+tE(km s -hV) 

12 v-^ cos(q x x) cos(q v y) cos(q z z) _ , , 

+A/ W — }J ^^Qq- (36) 

q q 

The system's wave function <£>(r, Qq) has to contain coordinates of an electron r and ions' displacement Qq. Usually 
they suppose the wave function is product of two parts depending only on electron coordinates and phonon coordinates 
accordingly: 

$(r,Qq) = *(r)0(Qq + <5Q q ) (37) 

*W=(^) 3/2cx p(-^)' (38) 

where the electron wave function is chosen in Gauss view to describe a localized electron in a cluster, ro is a polaron 
radius playing a role of a variational parameter. The phonon wave function <p is a wave function of a harmonic 
oscillator centered on the equilibrium displacement — <5Qq (an electron deforms a cluster and ions displace in new 
centers of equilibrium), which must be determined. Let us average the hamiltonian (|36p over electron's coordinates: 
H(Q q )= J^(r)H^{r)d 3 r: 

* «q) - T E W + l^ql 2 ) + 7^2 + E L ^ ( 39 ) 
q * o q 

where 

* - Mo\fi r r r vu™^™^™^^. m 

V L A J- L / 2 J- L / 2 J-L/2 <1 

The first term in the hamiltonian (|39[) describes harmonic oscillations (optical phonons) about equilibrium positions 
Qq = (zero oscillations and excited phohons if they are). The second term 3?i 2 /4m*rQ is a kinetic energy of an 
electron in the localized state. The last term is a potential energy of an electron's interaction with a deformation 
field. Following [l6| let us suppose the equilibrium displacement is £Qq = Lq/huj. Then the term in the electron's 
energy which is linear in Qq disappears. This gives us the energy of a system electron+dcformation field I(tq): 

* Wq) = ^ E (Wq + ^ql 2 + l^ql 2 ) + ^ + i E L \ (4D 



E(IQq + ^ql 2 + l p ql 2 ) +/ M- ( 42 ) 
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The first term in (|41l) describes the harmonic oscillations about new equilibrium positions — SQq. The second term is 
a kinetic energy of an electron in the localized state. After the done transition in the hamiltonian the last term 
is a potential energy of an electron's interaction with a deformation field plus an energy of the deformation field. 



/(eV) „. 




Figure 7: The graphics of a dependence of a system's energy electron+cluster (|17|) on a polaron radius I = I(ro) for an ionic 
cluster consisting of NaCL. The curves were build for edges of cubic clusters: 3a, 5a, 7a, 21a, 49a, where a = 2.81A is interatomic 
distances. The bold line is an energy of an infinity crystal with an electron. 
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Figure 8: The energy of a system electron+cluster (consisting of NaCl) as a function of a diameter of the cluster I(L) (solid 
line). The energy of a polaron in an infinity crystal L = oo is shown by a dotted line. Lengthes of edges (in interatomic 
distances a) of cubic clusters n = 3, n = 5, n = 7 are marked by vertical lines. 



Summation over q is equivalent to summation over n x ,n yi n z in accordance with Eq. (|30p . Since the summation is 
spreading on odd n only, then it is convenient to passing to new variables: n x = (n x — l)/2,n y — (n y — l)/2,n 2 = 
(n z — l)/2. Besides in the formula (|40[) the limits of integration can be expanded to ±oo that simplifies the calculation 
without a essential error because the restriction (|30[) on wave vectors of phonons does a main contribution in the effect 
of a finite volume of a cluster. Hence the energy can be written in a view: 



I(ro) 



3h 2 Ml 2 
4m*rn 2Huj L 3 



y,(n-l)/2 y. 



(n-l)/2 
=(— n 

(n-2)/2 



(n-l)/2 
ra 8 =(- 

(n-2)/2 



S^Z" for 
(-n+l)/2 ^n 8 =(-n+l)/2 ' 



odd 



— n/2 n y ——n/2 n z ——n/2 



yj(«-2)/2 

L 2 exp (-f^^ [(2S X + l) 2 + (2H y + l) 2 + (2n z + l) 2 ]) 



,for 



(2n, + l) 2 + {2n y + l) 2 + {2n z + If 



(43) 
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For an infinity crystal (n — > oo) the system's energy and the polaron radius are: 



i(ro) 



3h 2 Ml 1 



(44) 



mMo 



The results of calculation of /(ro) for different n are shown in Fig.tJTJ). We can see for clusters with a length of 
an edge n > 5 (n is number of interatomic distances a = 2.81A) a bound state of an electron exists, moreover the 
polaron radius nearly does not depend on a size of a cluster and it is equal to the polaron radius in an infinity crystal 
(]44fl tq m 2.8A. With increase of a cluster's size an energy of a polaron aspires to its energy in an infinity crystal 
/oo = — 0.26eV. Compared Fig. (|7) with the results in Fig.([T|) obtained in the quasiclassical approach we can see the 
quantization of cluster's oscillation brings to the same critical diameter of a cluster 5a. 

The bulky expression (|43p containing a triple sum can be simplified by a continuous approximation. The polaron 
radius ro nearly does not depend on a cluster's size L and it is equal to its value in an infinity crystal. Then in the 
expression (I43|) we can suppose ro equals to the optimal value (|44[) . Hence the energy of a cluster with an electron is a 
function of the cluster's size only: / = I(L). In a cube cluster with a length of an edge L a phonon with a wave vector 
q < ir/L cannot propagate. The upper limit of wave vectors' values caused by a crystal lattice does not influence on 
the result essentially. Then the cluster can be supposed spherical with a diameter L and we can replace the triple 
sum (|4*3"f by an integral: 



A result of the calculation of an energy of the system electron+cluster with help of the formula (|45|) is shown in 
Fig.©. Let us compare the result with an exact result shown in Fig.®. We can see the critical sizes of a cluster 
in both methods of calculation coincide: L C r = 4a -j- 5a, energies of the system for corresponding sizes of a cluster 
are nearly equal. At increase of a cluster's size to infinity the energy of a polaron aspires to its value in an infinity 
crystal. These facts afford ground to use the simplified expression (|4"5)) instead of the bulky formula (j4"3"f . 

IV. THE QUANTUM APPROACH. A WEAK COUPLING THEORY OF AN ELECTRON IN A 

CLUSTER. 

Localization of an electron in a cluster (crystal) because of an interaction with optical phonons can be described 
by a perturbation theory. This supposes a mass operator or an effective mass of an electron have peculiarities at the 
localization of the electron. Further application of the perturbation theory for a system with values of parameters 
(electron phonon coupling constant, size of cluster) exceeding the critical values is impossibly. 

Let Go and G are a free propagator and a dressed propagator of an electron accordingly (here % = 1): 




(45) 



G (k,e) 



1 



G(k,e) 



1 



(46) 



8 — 77- — h iS 1 

2m 



G - 1 -E(k, £ ) 



where 6 —> + , S(k, e) is a self-consistent mass operator: 




(47) 



A)(q, ui) is a propagator of a optical phonon: 



A)(q,w) 



2w 



(48) 



.2 



We suppose optical oscillations are dispersionless u 7^ u (q). After integration in (|4"Tf over the energetic parameter oj 
we have: 
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If in a crystal one electron propagates only then a corresponding quasi-particle (dressed electron) coincides with 
the particle (in a many-particle system quasi-particles are not identical to particles of the system). This means the 

2 

quasi-particle has an infinity lifetime, hence it propagates on a mass shell: e = + S(k, e). This equality brings to 
satisfaction to the self-consistence conditions (HHJ). 

For a crystal the summation over q is done over all vector space. As consequence of the boundary conditions 
(I27I29P on a surface of a cluster the summation must be done over the wave vectors (|30)l only. The system must 
be homogeneous so that a notion of the mass operator £ has a sense. But one cluster in any medium is a spatially 
inhomogeneous system. We propose a method shown in Fig.(|9]). Let us divide an ionic crystal into cubes with an 
edge L (a size of the cluster). Then let us introduce the boundary conditions (|2T|) on each cube faces as on a cluster's 
surface - a potential of the polarization field tp — 0. Let an electron propagates in this quasi-homogeneous system and 
interacts with optical phonons. The boundary conditions are imposed on the optical phonons. Then we can calculate 
the mass operator because a diagram series can be uncoupled (a momentum is conserved). 
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Figure 9: A division of the crystal into cubes with an edge L. In the system a charge e is propagating with a momentum k. 
On each cube face the boundary condition (|27|) has been given - a potential of the polarization field ip = 0. 

In a section (|ILI|) we showed a summation over the vectors ([50)) can be replaced by an integration over q from tt/L 
to +oo. Then the mass operator has a view: 

e = r dx r d q 



27rm k r 1 2mdx ( tt/L — kx \ 

= — arcsm —=^ — I arctan . (50) 

k \J2mujQ J i \/2mu)Q — k 2 x 2 \ \/2muj a — k 2 x 2 J 

The second term determines the finite size correction to the mass operator = — arcsin -^=== . Supposing 
L — > oo we have S = Sqo. 

A dispersion law of a quasiparticle is determined by a pole of the propagator G: e = + E(fc). For the limit 
k — > we have an expression: 

tt/L \ 









— OiUjQ 


H 


arctan 






k 2 










H 


+ 2m 







tt/L 2yJ2mujo tt 2 / L 2 — 2mwo 



TV 

where we introduced an effective mass of an electron m* : 



y/2muo L (tt 2 / L 2 + 2muioy 

2 tt/L \ k 2 

au; 1 arctan -==)+- — , (51) 



1 c I 7r dlCTdn ^JJE L (7r2/i 2 +2m Wo ) 2 



(52) 



From ([5T|) we can see a bottom of a conduction band falls. At L —> oo a fall of the bottom correspond to a result for 
an infinity crystal —au. The expression for an effective mass (|52p is more interesting. In an infinity crystal we have 
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the expression to* = to/ (l — |n. At a = 6 the effective mass turns into an infinity to* (a = 6) = oo, that means the 
autolocalization of an electron. At a > 6 the perturbation theory is not applicable. If the crystal is finite (cluster as 
in Fig|9|) then the effective mass is a function of a cluster's size m*(L). In an above-mentioned domain of values of a 
electron-phonon coupling constant such a value L C r exists that m*(L C r) = oo, to < m*(£ < Lcr) < oo, to*(0) = to. 
Thus if a cluster consists of a substance with a > 6 then a critical size of the cluster exists when an electron can be 
autolocalized. A dependence of a electron's effective mass on a size of a cluster consisting of NaCl is shown in Fig. ([T0| . 
The term ^ arctan gives a main contribution in the dependence (|52p because 

,. 2 7r/£ , 2 vr/L 

hm — arctan — ^=^= = 1, hm — arctan — ^=^= = 

L-yO 7T v / 2TOWo L-voo 7T v / 2toOJo 

2V2tow 7r 2 /i 2 - 2tow _ 2V2mw tt 2 /L 2 - 2muj _ 

L™ Z (7T 2 /L 2 + 2TOW ) 2 ' L^o Z (7T 2 /i 2 + 2TOW ) 2 ' 

that is just this term gives a required asymptotical dependence on L: 

£(£ — > oo) = Eqo, £(L^0) = 

to* (L — » oo) = to*(oo), to* (Z/ — > 0) = m (53) 

If to take into account this contribution only then we will obtain the same picture qualitatively with a not large 
divergence of the critical size. 




Figure 10: The dependence of an electron's effective mass on a size of a cluster (consisting of NaCl) m*(L)/m in a limit k — > 0. 
The critical cluster's size is Lcr ~ 39A ~ 14a. A graph on the insert is a dependence m*(L)/m caused by the main contribution 
— arctan S= only. In this case the critical cluster's size is Lcr ~ 28A ~ 10a. 

Above we considered a limit k — > 0. But the further consideration shows the effective mass is a function of an 
electron's momentum. In a domain k 2 /2m, < loq for an infinity crystal we have: 

k 2 „ .,. k 2 / a 3a k 2 15a 

\k) = -auj + — 1 



2m w 2m V 6 40 2mw 336 (2tow ) 2 



to*(/c — >• V2towo) = 7 v-. (54) 



VSarcsin^^-l 



1 - a (f 



In a domain k 2 /2m > ujq Cherenkov radiation of a phonon is possible. Let us generalize the expressions for an finite 
crystal with help of the expression (fSTJj) . A main contribution into the dependence on cluster's size L (a length of a 
cube edge or a diameter) is given by the expression ^ arctan ■ Then we obtain: 

EfX) = Eoo | 1 - - arctan ^2= ) , 

V 7T V2TOWQ / 



14 



m 



L(A) 



Figure 11: The dependence of electron's effective mass on a size of the cluster m*(L)/m in a limit k — > \f2muJo. The critical 
cluster's size is Lcr ~ 6.8A « 2a 4- 3a. 



m*(fc — > 0) = 



9 fc 2 

20 2muo 



m*(fc — > V2ma;o) 
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a- 1 



56 (2mw ) 1 



I- a § -1 



1 — — arctan 



tt/L 



— arctan Z/ L 

it \f Amu® 



(55) 



A dependence of an electron's effective mass on a cluster's size in a limit fc — > ^/2tocl>o is shown in Fig. pTj) . At these 
momentums of an electron a more strong coupling with optical phonons takes place. Hence the critical diameter of 
a cluster is much less than one in the limit k — > 0. Thus the weak coupling approach gives a dependence of the 
critical size on an electron's momentum Lcr(fc)- However the critical sizes in this approaches are several interatomic 
distances (oc 5a ~ 10a) as in the strong coupling approach. 



V. CONCLUSION. 



In this article we calculated a polaron's energy in an ionic cluster (by the example of a nanoparticle consisting of 
NaCl) and a critical size of the cluster regarding in formation of a autolocalized state of an additional electron. The 
calculation is done by both the quasiclassical method and the quantum approach (in a sense of quantization of a 
deformation field). Moreover we generalized a weak coupling method for a case of a polaron in a cluster. The results 
of the work are following. 

As in numerically solution and computer simulation in Q we obtained that the bound state of an electron in clusters 
from polar substances can be realized starting with a some critical radius in consequence with the boundary conditions 
on a cluster's surface for a potential of the polarization field <p(R) = 0. With increase of a cluster's size the binding 
energy increases and aspires to a polaron energy in an infinity crystal. A radius of an electron's autolocalization in 
a cluster (a polaron size) is such that an energy of a system electron+polarized cluster is minimal. However by the 
analytical variational method we found that this radius nearly does not depend on the cluster's radius. The binding 
energy increases monotonically with increase of the cluster's size. Moreover the pattern is observed: the bound state 
exists while the polaron radius is less than the cluster's radius. That is an autolocalized electron is hidden in a depth 
of the cluster and has a negligibly low probability to be on the outside. If the cluster is enough small so that an 
electron can be on the outside the bound state disappears. Unlike results of numerical simulation in 2, 5], Q where 
surface states were found which can be understood as excited state we found that if the polarization well is enough 
deep (cluster's size is enough big) then except s-state of the autolocalized electron other discrete levels are possible. 
We calculated energy of p-state and obtained a cluster's critical radius for this state. It is necessary to notice that 
similar results ware obtained in but another functional I(^,(p) was considered and ground state energy of an 
electron in a cluster and the autolocalization radius were expressed by static dielectric constant only. Hence in this 
model autolocalisation effect is not polaron effect unlike our model. 

An interaction of an electron with the polarization field is the interaction with longitudinal optical phonons. Due 
an ionic type connection between atoms ionic clusters form stable cuboid structures and have a cubic microlattice. 
However it is necessary to notice for clusters with a size bigger the the critical size R 3> fo an excess electron has 
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low probability to be near a surface of the cluster hence the shape have not principal importance. The boundary 
condition on a cluster's surface for the polarization field <p(R) — brings to a restriction on possible wave vectors of 
optical phonons n/L < q < oo. The optical phonons are standing waves between opposite sides of a nanoparticle. 
We obtained the polaron energy in a cluster (|43[) taking into account the restrictions on a propagation of phonons. 
Results for an energy and a radius of the polaron correspond to the results of the quasiclussical approach - existence 
of a critical size of a cluster regarding in formation of an autolocalized state of an electron, an asymptotical aspiration 
of a binding energy of an electron in a cluster to the binding energy in an infinity crystal with an increase of a 
cluster's radius, an independent of a polaron radius on a cluster's size, an infinity low probability for an electron to 
be on the outside of a cluster if the electron is in the autolocalized state. These results are different from the results 
0-Q of computer simulation of electron localization in alkali-halide clusters M+X~_ l and M+AT~_ 2 which contain 
F-center defects and an electron localized near halide ion. Unlike our result, where an excess electron autolocalizes on 
a neutral ionic cluster due interaction with the induced polarization, such systems have not any critical parameters. 
This problem is analogous to a problem of F-center model of negatively charged metal-ammonium clusters [3] . It is 
necessary to notice that an essential difference of a polaron state in a cluster from the polaron state in a crystal is 
localization of an electron in the center of a cluster. On the contrary in an infinity crystal all points are equivalent 
(in a continual representation of a medium) hence the electron can move. 

We described an electron's localization in a cluster due an interaction with longitudinal optical phonons by a 
perturbation method taking into account the boundary conditions for the electric field of a polarization. The electron 
localization in a cluster manifests itself as a singularity of an electron's effective mass m(L CT ) = oo. Using the 
perturbation method we showed a dependence of the critical size of a cluster on a electron's momentum in a conduction 
band L cr = L cr (k). The largest critical size of a cluster is for an electron situated on a bottom of the conduction 
band, the smallest critical size is for an electron's momentum on a limit of Cherenkov radiation of an optical phonon. 
This fact is explained by the fact that an interaction of an electron with an ionic crystal increases with an increase 
of an electron's kinetic energy, because the larger electron's energy the stronger an electron radiates virtual phonons 
and then interacts with them. 



[1] P. Stampfli, Bennemann K.H. Phys. Rev. Lett. 58 (1987), 2635 

[2] V.D. Lakhno , Chem. Phys. Lett. 240 (1995), 585; N.K.Balabaev, V.D.Lakhno, A.M.Molchanov, B.P.Atanasov., J. Mole- 
qular Electronics 6 (1990), 155; N.K.Balabaev, V.D.Lakhno., Zh. Fiz. Khimii 65 (1995), No8, 1358; V.D.Lakhno, A.N. 
Pancratov, Bulletin of the Russian Academy of Sciences: Physics 64 (2000), No8, 1172. 

[3] V.D.Lakhno, Bulletin of the Russian Academy of Sciences: Physics 62 (1998), No6, 885. 

[4] R.N. Barnett, U. Landmann, G. Markov, A.Nitzan, J. Chem. Phys. 93 (1990), 6226 

[5] R.N. Barnett, U. Landmann, A.Nitzan, Phys. Rev. Lett. 62 (1988), 106 

[6] R.N. Barnett, U. Landmann, C.L. Cleveland, Phys. Rev. Lett. 59 (1987), 811 

[7] U. Landman, D. Scharf, J. Jortner, Phys. Rev. Lett. 54 (1985), 1860 

[8] E.C. Honea, M.J. Homer, P. Labastie, R. Whetten, Phys. Rev. Lett. 63 (1989), 394 

[9] G. Rajagopal, R.N. Barnett, A. Nitzan, U. Landman, E.C. Honea, P. Labastie, M.J. Homer, R. Whetten, Phys. Rev. Lett. 
64 (1990), 2933 

[10] SeGi Yu, V. B. Pevzner, K. W. Kim, and Michael A. Stroscio Phys. Rev.B 58 (1998), 3580 
[11] J. J. Licari, R. Evrard, Phys.Rev.B 15 (1977), 2254 
[12] Stroscio M.A., Phys.Rev.B 40 (1989), 6428 

[13] de la Cruz R. M., Teitsworth S. W., Stroscio M. A., Superlattices and Microstructures 13 (1993), 481 
[14] D.O. Welch, O.W. Lazareth, G.J. Dienes, R.D. Hatcher , J. Chem. Phys. 68 (1978), 2159 
[15] A.S. Davydov, Theory of Solids, Nauka, Moscow (1976), p. 640 (in Russian). 
[16] Gerald D. Mahan, Many-particle physics, Third edition, p. 785 



